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BURGERS APPROXIMATION FOR TWO-DHIENSIONAL 
FLOW PAST AN ELLIPSE 


By 

J. M. Dorrepaal* 

SUMMARY 

This paper examines a linearization of the Navier-Stokes equation due 
to Burgers in which vorticity is transported by the velocity field corre- 
sponding to continuous potential flow. The governing equations are solved 
exactly for the two dimensional steady flow past an ellipse of arbitrary 
aspect ratio. The requirement of no slip along the surface of the ellipse 
results in an infinite algebraic system of linear equations for coefficients 
appearing in the solution. The system is truncated at a point which gives 
reliable results for Reynolds numbers R in the range 0 < R < 5. 

t 

Predictions of the Burgers approximation regarding separation, drag and 
boundary layer behavior are investigated. In particular, Burgers lineariza- 
tion gives drag coefficients which are closer to observed experimental 
values than those obtained from Oseen's approximation. In the special case 
of flow past a circular cylinder. Burgers approximation predicts a boundary 
layer whose thickness is roughly proportional to This is in agree- 

ment with the nonlinear theory despite the fact that the Burgers calcula- 
tions are carried out using only moderate values of the Reynolds number. In 
the matter of separation, it is shown that standing eddies form on the down- 
stream side of a circular cylinder at R “ 1.12. Interestingly enough, this 
is the same value predicted by Skinner (1975) using singular perturbation 
techniques on the full nonlinear problem (see Van Eyke, 1975). 

The linearizations due to Oseen and Burgers both give spatially uniform 
approximations to the flow past a finite obstacle. The main difference is 
that vorticity is transported around the obstacle in Burgers flow rather 
than through it. The results of this paper verify the superiority of 
Burgers approximation in modelling the flow near the obstacle at low to 
moderate Reynolds numbers. 

*Assoclate Professor, Department of Mathematical Sciences, Old Dominion 
University, Norfolk, Virginia 23508. 
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1. INTRODUCTION 

» 

Owing to the formidable nature of the Navier-Stokes equation, the 
history of fluid mechanics research is filled with simplifying approxima- 
tions to this nonlinear problemt Included among these is a class of approx- 
imations which replaces the nonlinear inertial term (v • V)v by a linear 

one (v • V )v where v is given. Of particular interest are the Stokes 
0 o ^ 

and Oseen approximations where is constant. These have contributed 

significantly to the understanding of basic fluid dynamics behavior, espe- 
cially in the low Reynolds number regime where they are related asymptot- 
ically to the Naviut-Stokes solution. 
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A linearization which has received conaiderably leas attention ia 
Burger's approximation (1928) in which is taken to be the continuous 

potential flow around the body* As noted by Dryden et al. (1956), the 
vorticity equation in this csae is identical to the temperature equation 
used by Bousainesq in his study of the conduction of heat from a hot body 
placed in an irrotational fluid. In Burgers approximation the convective 

velocity field follows the surface of the body in its immediate 

neighborhood and approaches the velocity of a uniform stream at a great 
distance from the body. Burgers flow is asymptotically equivalent to Oseen 
flow far from the body, but in its immediate neighborhood Burgers flow 
models the exact flow more accurately. Of course Burgers approximation 

suffers from the defect that the convective velocity vector does not 

tend to 3 as one approaches the surface of the obstacle. Furthermore 
at moderate values of the Reynolds number R, separation occurs on the down- 
stream side of a bluff obstacle and the resultant velocity field v no 
longer resembles This limits the suitability of Burgers approximation 

to small values of R. Nevertheless Burgers linearization provides a 
spatially unifontv approximation to the solution of the Navier-Stokes equa- 
tion and it is the purpose of this paper to investigate the extent to which 
it improves upon the Oseen approximation. Attention is restricted to two 
dimensional steady flows past circles, ellipses and flat plates. 

2 . MATHEMATICAL FORMULATION 

The nondimens ional Navier-Stokes equation has the form 

R(v • ^)v =» -^p + V (2.1) 

where R, v, p are Reynolds number, fluid- velocity and pressure, 
respectively. The velocity must also satisfy the continuity equation 

^ • V “ 0 (2.2) 

which is guaranteed by the introduction of a stream function '^(x,y) 
defined by 


5 


(2.3) 


V " curl 

The problem for the stream function corresponding to (2.1) is given by 

[72 + RD(ij;)]aj - 0 (2.4) 


■ — U 


(2.5) 


where w(x,y) is the magnitude of the vorticity vector and 

D(t) . 

8x 9y 3y 9x 


Consider a finite obstacle whose boundary is denoted by B, Let 
be the stream function for the continuous potential flow around B which 
approaches a uniform stream at infinity. The boundary value problem for 
’f'o(*iy) is given by 


0 


( 2 . 6 ) 


'I'O 


0, ipo ~ y as x2 + y2 


(2.7) 


Burgers linearization is defined in the following way; 
[y2 + RD(\|)o)]w ■ 0 


( 2 . 8 ) 


72^ a -(j) 

♦I -I! 


0 , ~ y as x2 + y^ -► 


(2.9) 

( 2 . 10 ) 


Equation (2.8) is called the vorticity equation and (2.9) is Poisson’s equa- 
tion. The Oseen linearization is obtained by substituting i|>o(x,y) » y in 
( 2 . 8 ). 
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3. SOLUTION OF THE VORTtCTTY EQUATION 


The exponential decay of the vorticlty in Oseen flow suggests that we 
attempt a solution of (2. 8) of ':he form 


w(x,y) « F(x,y) exp [f^(x,y)] 


(3.1) 


The substitution of (3.1) into (2.8) gives a second order partial differen- 
tial equation for F(x,y) with non-constant coefficients involving deriva- 

3F 3F 

tives of i|) and f . If we set the coefficients of — and — to zero, we 
obtain the following: 


3f^ , S'I'O 


3 X 


2 ^ 3y ' 


3f 


3y 


3ij^ 

2 3x 


(3.2) 


These are the Cauchy-Riemann equations. They suggest that we choose 
f^(x,y) to be i R times the velocity potential of the irrotational flow 

past B. The resulting equation for F(x, y) is given by 


F - — q(x,y)^ F » 0 

4 


(3.3) 


where q(x,y) =■ 


3ijj 


o 


’3ij; 


+ 


1/2 


is the magnitude of the potential flow 


k.3 X / \9 y / 

velocity. In Oseen flow we have q(x,y) = 1. 

Consider a curvilinear coordinate system (5 ,h ) defined by 


X + iy = f (? + ih ) 


(3. 4) 


where f is an entire function. The metric coefficient a(C,n) for this 
transformation is defined by 




h ) W 


3x\2 ^/dyV 


3n 


3n 


(3.5) 
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The vorticity in (C >n ) 
w(5,n) “F(?,n) 


coordinates is expressible in the form 

exp [i R|i (5 ,r \ )] 

2 


(3. 6) 


where is the harmonic conjugate of \(>o. From (3* 3) the equation for 
F(5 ,n ) is given by 


9^F ^ 3^F 
3?^ 3n^ 



p (5 ,n P - 0 


(3.7) 


where p (I ,n ) 





1/2 


a(5)n) q(5>n). Equation (3.7) is 


L\95 / \3n / J 

separable provided P^ is exnressible in the form 

■ P(5,n)^ » g(5) + hOf) (3:8) 

where g and h are arbitrary functions. Under this assumption we have 
F(?,n) » E(5) H(n) (3.9) 


where 


E” - [o + i g(5 )] E » 0 
4 


H" + [a - i h(n )] H = 0 
4 


(3.10) 


(3.11) 


and a is the separation eigenvalue. 

There are three geometries in which the vorticity equation can be 
solved exactly. Consider first the case of a circular cylinder of unit 
radius. The appropriate function in (3. 4) is 

f(z) “ e^ (3.12) 
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which implies that x •• cos hi y • e^ sin n. The connection with 
polar coordinates Iv r “ e^, 0 ■" n and the unit circle is given by 5 » 0, 
The potential flow boundary value problem is given by? 





_o 

3n2 


n 0 


o,u) 


toC0»Tl) " 0» ~ n as 5 *»• »* 


The solution is 

“ 2 sinh % sin n 

^(C>n) ■ 2 cosh 5 cos n (3.14) 

from which we obtain 

p(5,n)2 » 2 cosh 25-2 cos 2n (3.15) 

Thus the separability condition (3.8) is satisfied and the corresponding 
separated equations are, from (3.10) and (3.11), 

E" - [a + - r 2 cosh 25 ] E - 0 (3.16) 

2 

H" + [a + i r2 cos 2ti] H » 0 (3.17) 

2 

These are respect- ively the modified and the conventional Mathieu equations. 

A brief discussiosj of their solutions is included in Appendix I. 

In a streaming two dimensional flow past a symmetric body, the vortici- 
ty must be odd in h and periodic with period 2iT. It must decay exponen- 
tially as 5 00 , except possibly in the wake n ■ 0, The periodicity condi- 
tion determines the eigenvalues in (3.17). The corresponding odd 

eigenfunctions are denoted by se (ri, - — R^} . In (3.16) the eigenfunctions 

4 
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which decay exponentially are denoted by Gek (5 , “ il R^l • 11108 from (3.6), 

^ 4 

(3.9) and (3.14) the vorticity function has the form 

00 

w(5 ,n ) - - exp [R co8h 5 cos n ] I W GeR C? > - i R^ ) 

n-1 " 4 

• sefn, (3.18) 

4 

where the minus sign is included for convenience and the coefficients 
are constants to be determined. The aaymptotic behavior of the vorticity 
is, from Appendix I, 

w (5 ,h ) ~ w(n ) exp [ ~ i ? -- i R e^ (1 cos r| )] 

2 2 

» w(0 ) r”^/2 exp [ - i. Rr (1 - cos 6 )] (3.19) 

0 

where w(0) ■ 0, w'(0) ^ 0 and w(9 ) has period 2ir. For purposes of compari- 
son we have solved the vorticity equation for Oseen flow past a circle in 
Appendix II and it is observed that the Oseen vorticity has the same 
asymptotic behavior. 

Consider next the flow past an ellipse where the major axis of the 
ellipse is parallel to the flow at infinity. The appropriate coordinate 
system (? ,n ) is defined by 

X + iy ■ a cosh (5 + in ) « a cosh 5 cos n + i a sinh 5 sin n (3*20) 

t 

The curve 5 "Co is an ellipse with major axis of length 2a cosh Kq and 
minor axis 2 a sinh in order that the unit of length be the semi-major 

axis, we choose 

a ■ sech ? (3.21) 

0 

It is convenient to make the transformation 

? "5-5 (3.22) 

o 
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Thus from (3.20), the coordinate system (C,n) is defined by 


X + iy » a cosh (C + 5 + in) 

o 


and the ellipse is given by C “ 0. The metric coefficient for this trans- 
formation is a(?,n) where 

0(2 ■ i a^ [cosh 2(C + 5^) - cos 2n] (3.: 


The potential flow stream function '|'o(?>n) satisfies the boundary value 
problem 

V2 i|, » a“2 °+ 2. -0 (3.2! 

° 3^2 3^2 

1 ^o"*^ ^ . C 

'<jf (0,n) = 0, ^ (?,n) ~ ij s e sin n “ A sin n as C “ 


1 Cq 

whe,re A =» — a e . The solution is 


il>^(5,n) “ 2A sinh 5 sin n 
<}>(?, n) =* 2A cosh C cos n 


Because of the similarity between (3.14) and (3.26), the solution of 
the vorticity equation in the elliptical case follows the method already 
outlined. Thus the vorticity function has the form 


w(c,n) “ - exp [AR cosh ? cos n] I W Gek (?, 

n»l ^ " 

• se^[n, - — a2r2) 

4 


i A2 r2) 
4 


The Reynolds number R is based on semi-major axis as length unit. The 
circular cylinder result in (3.18) is recovered by letting * and 

observing that 
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(3.28) 


A 



sech 


'o * 1 


If the ellipse is oriented so that its major axis is perpendicular to 
the flow at infinity, we use the coordinate system (X,n) defined by 

X + iy " a ainh (X -f in) (3.29) 

The ellipse is given by X " Xq. As before we define a modified coordinate 
system (y,n) by letting p " X - Xq. The metric coefficient is a(p,n) 
where 

jj2 2(y + X ) + cos 2n] (3.30) 

A O 


The potential flow stream function '|’o(h>n) satisfies a boundary value 
problem identical to that given in (3.25) and so the solution of the 
vorticity equation proceeds exactly as before. The vorticity for this case 
is obtained from (3.27) by making the following adjustments: 

C •<- p 

1 ^ 

A « — e ° sech X (3.31) 

ft o 


4. SOLUTION OF POISSON'S EQUATION 


Poisson's equation (2,9) can be solved using a Green's function 
approach. Equation (3,27) gives the form of the Vorticity function for the 
three geometries considered in section 3 and so the problem reduces to the 
following; 


= l«(C,Ti)r2 


3^ 

95^ 



w(C,n) 


'P(o,n) 


0 , 


i|^(C,n) A e*" 


sin n as 5 


( 4 , 1 ) 
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The appropriate Green’ s function is defined by 

G(c ,ti Ic ’ ,n ' ) - 5(n - n’) 

G(0,n k ' ,Ti ' ) o» G(<» ,n k 'h ’ ) - 0(1) 

The solution of (4. 2) is given by 

Oft ,n k ■ ,n ' ) - - ^ ^ -..Jill 

4u Uosh (? + C ' ) - cos (n - n ' ) J 


(4.2) 


(4.3) 


Proceeding in the usual way, we have 

'|'(C'>n') - // [CSV^ij; - ijiV^G] d?dn + // wGa^ dCdn (4.4) 

A A 

where A is the fluid region between the obstacle ? « 0 and the curve C « P 
with r being large. 

Green's Identity transforms the first integral in (4.4) into contour 
integrals around the boundary curves. The Integral around ? “ 0 vanishes 
because of the boundary conditions on G and Thus we have 


// [G72t|^ - i|>v2g] d?dn 

A 

= /’ ( 
"IT \ 

0 - 1 ii'l dn 

^ 'Vc = r 

(4.5) 

Asymptotic expressions for 

G,i|^ 

and their derivatives 

are given below 

~ A e^ sin n + 0(1) 



(4.6) 

9tlJ C . 

— ~ A e sin n + 0(e ) 

3? 



(4.7) 

G ~ i- c ' + i sinh ? ' e~^ 
2ir 7T 

cos(ri 

- n ' ) + 0(e ^^) 

(4. 8) 

a G 1 . . , , ■% 

— ~ - _ sinh ? ' e cos 

(n “ n 

’ ) + 0(e ^ ) 

(4.9) 


9? 7T 
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Substituting these into (4.5), we obtain the following result; 


lira g(^ - '1'-^) dn “ 2A sinh 5 ' sin n‘ ■ (4.10) 

r-Ho -7T \ac - r ° 


Thus in (4.4) we have 


iKc',n') - ’!' (?', n*) + r /” a^(c,n) <jj(c,n) G(c,nl?',n')dcdri (4.11) 

"u o 


5. DETERMINATION OF VORTICITY COEFFICIENTS 

The only boundary condition from (2.10) which remains to be satisfied 
is the no-slip condition. Its invocation yields unique values for the 
vorticity coefficients (W^). Thus in (4.11) we require 

(0,n') » 0 (5.1) 

H' 


which implies that 


a^o) d;dn (5.2). 

“tt 0 cosh C - cos (n - n') 


(1965), we have 


“1 + 2 'I e cos k(n - n') (5.3) 

cosh Z ■* cos (n - n') k“l 

But a^(?,n) w(5,n) is odd in n and so the even part of (5.3) makes no 
contribution to (5.2). After substituting (5.3) into (5.2) and equating 
coefficients of sin kn' to zero (k * 1,2,3, ...), we have 

f“ w e sin kn d?dn ®* -2 Att5j^^ k = 1,2,3,... (5.4) 

"ir o 


0 “ 2A sin n' + — 
2ir 

From Gradshteyn and Ryzhik 
sinh C 
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The substitution of (3. 27) into (5. 4) yields the following infinite linear 
system for the unknowns {Wj^} 


I Wn^nk"^'"^kl 1,2,3,... (5.5) 

n"l 


where 


7T ^ 

^nk " ^ cosh C cos n 1 “(? >n)^ 

o 0 

• Gek (? , - i ) se (n , - i ) sin kn d?dn (5. 6) 

"4 ”4 

To obtain an approximate solution to (5.5), we truncated at n * k ■ 8. 

This necessitated the calculation of 64 integrals of the form (5. 6) which is 

the chief drawback of this procedure. Nevertheless we were able to obtain 

good results in the range 0 < R < 5. Tn calculating the coefficients 

C it is convenient to choose a modulus q for which the eigenvalues 
n cC ^ 

On in Mathieu’s equation are tabulated. The corresponding Reynolds 
number can be found from the relation 


R « 2A"^ |q 1^^^ (5. 7) 

We examined five geometries: the circle, the ellipse with aspect ratio 3:1 

oriented parallel and perpendicular to the flow and the flat plate oriented 
parallel and perpendicular to the flow. The first eight vorticity coef- 
ficients for a sampling of geometries and Reynolds numbers are given in 
Tables 1-5. 

The rapid decay of the coefficients W^ does not mean that the 
vorticity series (3.27) converges rapidly in all regions of the flow. When 
q “ -1, for example, Gek (0,-1) = 0(106) and Gek (0,-1) = 0(1010) and 
the decay of the W^^’s is offset by the growth of the function values 

(0>”1)* The vorticity series may converge rather slowly therefore in 
the vicinity of the obstacle and indeed this is observed in the elliptical 
cases. The rate of convergence appears to be a maximum in the case of the 
circle and we therefore investigated separation phenomena for this geometry 
only. 


15 


Table 1. Circle. 


Table 2. 



R - 2 

.0 

R “ 4.0 



13. 250 


773. 03 


Wa 

-2. 2884 


-617.12 


W3 

0. 37690 X 

10-1 

11.281 


\ 

- 0 . 24210 X 

10 "^ 

-0.71475 


W5 

0. 62517 X 

10-6 

0.12619 X 

r—4 

1 

0 

\ 

-0. 78300 X 

10 -^ 

-0.11800 X 

10 -^ 

W7 

0. 53897 X 

10-12 

0. 60380 X 

10-6 

Ws 

-0.21073 X 

10 -^® 

-0.15441 X 

!-■ 

0 

1 

03 


Ellipse (aspect ratio 3:1) parallel to flow = i An 2] 

^ o Z ' 



R = : 

i. 

12 

R = 4.24 

% 

5.1649 



103. 41 



-0.41325 



-46. 711 


W3 

0. 62398 

X 

10-2 

1. 8916 


\ 

-0.46197 

X 

10 -^ 

-0.64596 

10-1 

W5 

0. 20407 

X 

10-6 

0.11379 X 

10-2 

Wg 

-0.56252 

X 

10 “® 

-0.12250 X 

10 -*^ 

W 7 

0.10014 

X 

10-11 

0. 78668 X 

10-7 

Wb 

-0.10497 

X 

10 -^*^ 

-0. 25930 X 

10 “® 


X 

Table 3. Ellipse (aspect ratio 3:1) perpendicular to flow (x “ 4 An 2) 

O M 



R - 2 . 

12 

R - 4.24 

Wi 

5. 8811 


82.188 


\iz 

-0.33473 


-30. 625 


^3 

-0.10127 X 

10”2 

0. 13610 



0.25251 X 

10“^ 

0.15274 X 

10“^ 

^5 

0. 27597 X 

10-7 

-0.11094 X 

10-3 

We 

-0.30245 X 

0 

1 — 1 

-0. 23548 X 

10-^ 

W 7 

-a 13758 X 

10-12 

0.14100 X 

10-7 

We 

0.83636 X 

10-15 

0, 22353 X 

0 

1 

1 — • 

0 


Table 4. Flat plate parallel to flow “ 0) 



R » : 

1, 

83 

R - 4.0 


Wl 

6. 7127 



25. 215 



-0.58416 



-5. 2619 


'^3 

0.11229 

X 

10-1 

0.17616 



-0.10123 

X 

10-^ 

-0. 32847 X 

10-^ 

W5 

0. 55440 

X 

10-6 

0. 34490 X 

10-^ 

% 

-0.17737 

X 

10-® 

-0. 22002 X 

10-® 

W7 

0. 35771 

X 

10-11 

0. 82930 X 

10-3 

Ws 

-0.41796 

X 

3^0^ 4 

-0.17796 X 

10-^ ^ 
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Table 5. Flat plate perpendicular to flow (Xq ■ 0) 



R » 2. 

83 

R - 4.0 


Wi 

9.0389 


27. 223 



-0. 40074 


-3. 0988 


^3 

-0. 73015 X 

10-2 

-0. 74847 X 

10-1 


0. 62359 X 

10-*+ 

0.15255 X 

o 

I — i 


0. 34806 X 

10-6 

0.13029 X 

10-^ 

We 

-0.13098 X 

10-^ 

-0.12114 X 

10-^ 

Wy 

-0. 25803 X 

10-n 

-0. 37198 X 

10-9 

Wb 

0.62502 X 

10“^ 

0.21338 X 

10-n 


6. SEPARATION 

Separation occurs on the downstream side of the circular cylinder 
provided the l^ynolds number exceeds a critical value defined as the 

value of R for which 

H (0,0) - 0 (6.1) 

3n 

For Burgers flow we find R^ * 1.12 which is a new result. Yamada (1954) 
has shown that R^, » 1. 51 for Oseen flow and Underwood (1969) has obtained 
=* 2. 88 from a numerical solution of the full nonlinear equation. 

The Burgers result should be less than the numerical value. The con- 
vective velocity field in Burgers flow is potential flow past the cylinder 
and this violates the no-slip condition at the cylinder's surface. The 
velocity field which solves the full Navier-Stokes equation satisfies this 
condition. Thus convection effects near the cylinder are more dominant in 
Burgers flow than in Navier-Stokes flow and any phenomena related to convec- 
tion, such as separation, should occur aL lower Reynolds numbers in Burgers 
flow. 
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The fact that the Burgers result is less than the Oseen value also can 
be explained. Separation begins at the rear stagnation point P of the 
cylinder where locally the flow appears as in Figure 1. At the onset of 
separation two eddies of circulating fluid form about P. (We refer to this 
pair of eddies as a separation vortex). The direction of motion along the 
axis of symmetry inside the vortex is opposed to that outside (Figure 2). 

In Oseen flow the convective velocity field is constant in magnitude and 
perpendicular to the cylinder boundary in the vicinity of P as shown in 
Figure 3. Oseen convection therefore deters the establishment of reverse 
flow at P because it directly opposes the direction of fluid motion along 
the axis of symmetry inside the vortex. In contrast the convective velocity 
field in Burgers flow vanishes at the point P and is small in magnitude 
near P (Figure 4). Burgers convection does not oppose the establishment 
of a vortex about P to the same degree that Oseen convec'.ion does and 
separation therefore initiates in Burgers flow at a lower Reynolds number. 

As R increases beyond R^ -1.12 in Burgers flow, the separation 
vortex grows in size, l^hen R » 2. 0 the flow appears as in Figure 5. The 
length of the vortex is PQ » 0. 53 where OP is the unit of length, and <S0P 
“ 34.8“. Point Tj whose 6 -coordinate is 83®, marks the location where the 
fluid pressure along the boundary is a minimum. The flow from T to S is 
against an adverse pressure gradient. 

7. CALCULATION OF DRAG COEFFICIENTS 


Several authors (Imai, 1951; Kawaguti, 1953; Dennis and Danwoody, 1966) 
have shown that an obstacle' s drag coefficient can be obtained from the term 
of 0(1) in the asymptotic expansion of the stream function. This obviates 
the calculation of stress components which often are difficult to obtain 
accurately. The task therefore is to find the leading terms in the asymp- 
totic expansion of (4.11). 

We first obtain a series representation for the Green's function in 
(4.3). From Magnus et al. (1966), we have 
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_ 0*\i^^-'^!* I *■ IW 

OF POOR QUALITY 


00 

I i e cos kx “ i t - i An (2 cosh t - 2 cos x) t > 0 
k-1 k 2 2 


Manipulation of (7.1) yields the following results; 


G(c ,n |5 ' ,n ' ) - i- c + i i i e"^*^ ' 

2tt tt k""! k 

• sinh k? cos k(n - n') C < C' 


I ^ ^ I r I -kC 
— ?'+-^ -e ’ 

2tt it k“l k 


• sinh k5 ' cos k(Ti - ri ' ) 5 > 5' 


Substituting these into (4.11) and simplifying, we obtain 


’1^(5 ',n') =« 4'^(?',n') + i I i sin kn' 

X k“l k 

• 0) sinh k? sin kn dCdn + — ~ sinh kC ' 

-IT o TT k"! k 

• /^ /* w e sin kn d?dn 

-IT C • 


An expression for the second integral can be obtained from (5.4): 


-7f -00 ^ 

J J w e sin kn dCdn = -2 Ait6 , 
-TT C' 

- /^ /^ u e sin kn dcdn 

-TT O 


Thus a series representation for the stream function is, from (7.4) 


t (C ' >n ' ) = i ^ i- e ^ s in kn ' /^ /^ to s in kn dC dn 
IT k=l k o o 


— I — .$in kn ' ct^ w c sin kn d?dn 

IT k=l k 0 o 


(7.1) 


(7.2) 
» 

(7.3) 

sin kn ’ 

(7.4) 


(7.5) 


(7.6) 
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This expression can be written in the form 


c * I 

'K?'»n') ■ A sin n' - — sin n' 


to e ^ sin h dijdri + trA 


o o 


-c 


- — 2 ~ sin kn' 

TT k-2 k 


w 

+ — ~ sin kn* 

IT k“l k 


to e sin kii.jd?dri 

o o 


e-^^' 

r / , k? 

/ / cr (ii Q sin kn dcdn 
o o 


M' 


(7.7) 


It can be, seen from (5.4) that the first two bracketed expressions in 
(7.7) are indeterminate forms ( ~ ] as «. The third expression is of 

the form |~| in the limit. By invoking A'Hopital's rule we obtain con- 
tributions from each expression to the 0(1) term in the asymptotic expan- 
sion of The result is 


K C ' » n ' ) '^ A e^ s in n ' + “ I ^ sin kn ' 

TT k“l k^ 


{lim to sin kn dn} + 0(e ^ ) 

5+00 o 


(7.8) 


The limit can be computed and the series summed. The details are included 
in Appendix III. The expansion has the form 


K 5' »0' ) ~ A e^ sin n' - — C (±1 - -2-) + 0(e ^ ) 

r» ^ 


(7.9) 


where, in the second term, the plus sign is chosen when 0 < n* < ir and the 
minus sign when -ir < n' < 0. This terra is analytic along n' “ ''f, but suf- 
fers a finite jump discontinuity along n* “ 0 which coincides with the wake. 
Dennis and Dunwoody (1966) comment that this terra must be present in order 
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to give non-zero drag. Kavaguti (1953) shows that the constant Cd is the 
drag coefficient for the obstacle in the flow. 

Table 6 contains drag coefficients for the circular cylinder in Os sen 

t 

flow. In table 7 the Burgers flow results are given. 


Table 6. Circular cylinder in Oseen flow. 



0 

1. 51 

2.0 

3.0 

4.0 

5. 0 

8. 

08* 

6.62 

5.85 

4.98 

4.50 

4.18 


D 


*From Toraotika and Aoi (1951). 

Table 7. Circular cylinder in Burgers flow. 




1.0 

1.12 

1.25 

1.50 

2.0 

2. 83 

3.46 

4.0 

7.76 

7.30 

6.86 

6.22 

5.34 

4.47 

4.04 

3.75 


Figure 6 is a plot of the results contained in Tables 6 and 7 along 
with some of Tritton’s (1959) experimental values. The graph indicates that 
Burgers approximation is an improvement over the Oseen approximation in 
modeling the flow past a circular cylinder. 

In Tables 8-11, a sampling of Burgers drag coefficients is given for 
elliptic geometries. The results are plotted in figure 7* 


Table 8. Ellipse (aspect ratio 3;1) parallel to flow. 


2.12 

3.0 

4. 24 

5.20 

3.93 

3.22 

2.65 

2.36 
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Table 9. Ellipse (aspect ratio 3:1) perpendicular to flow. 


R 2,12 

3.0 

4.24 

5.20 

5.35 

4.56 

3.91 

3.58 

D 





Table 10, Flat plate parallel to flow. 


R 

2.83 

4.0 

5,66 

6,93 

c. 

2.42 

1.94 

1.56 

1,34 


u 


Table ll, flat plate perpendicular to flow. 


R 2.83 

4.0 

5.66 

6,93 

C 4.54 

3.98 

3.54 

3.31 

D 





8, BOUNDARY LAYER THICK^?SSS 

The procedure outlined in section 7 can be extended to obtain higher 
order terms in the expansion of 'I'. When expressed in polar coordinates 
(r,6), this expansion has the form 

OP 

~ r sin 9 - i C (±1 - i) + I r"" sin n9 + x(r,9) (8.1) 

2 ° TT n»l " 

where the$j^ (n = 1,2,3, ...) are constants andx(r,9) is exponentially 
small in r. 
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Except along the line 0 ■ 0, the algebraic part of the expansion is harmonic 
and constitutes the potential flow far from the obstacle* 

Because of the tedious nature of the calculation, the expansions were 
computed for the circular case only* Both Oseen and Burgers expansions were 
calculated for purposes of comparison. The results are presented in Tables 
12 and 13* 


Table 12* Oseen flow. 


R ■ 

1* 51 

\j; ~ r sin 6 - 3.3118 ( 

±1 - i) - 0*5995 





R - 

2*0 

2. 9247 

0, 5777 

R " 

3*0 

2*4888 

0. 5485 

R - 

4. 0 

2* 2481 

0,5302 

R - 

5.0 

2. 0918 

0* 5171 




Table 13* 

Burgers flow. 

R “ 

1*12 

t ~ r sin 6 - 3* 6509 

(ii - i) 

— 1.1251 




\ ir / 


R " 

2* 0 

• 2* 6709 


1*1396 

R « 

2* 83 

2.2327 


1.1439 

R ■ 

3.46 

2. 0213 


1*14,86 

R * 

4.0 

1. 8763 


1*1467 


sin 0 
• r 


sin 9 
r 
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Since both Oseen and Burgers flows are spatially uniform approximations 
to the exact solution of the Navier-Stokes equation, they predict a boundary 
layer surrounding the obstacle. The outer edge of the boundary layer is 
defined to be the curve along which the algebraic part of the asymptotic 
expansion of vanishes. The curve so defined determines the displacement 
body which the potential flow far from the cylinder "sees-/" The displace- 
ment body includes the cylinder, its separation vortex and the surrounding 
boundary layer. 

By setting the expansions giveq in Tables 12 and 13 equal to zero, we 
obtain approximations to the displacement bodies for the various flows. A 
typical example is given in Figure 8. The displacement body is semi-infi- 
nite with its width at infinity numerically equal to the drag coefficient. 
Since the boundary of the circular cylinder is given by r » 1, the thickness 
of the boundary layer is easily calculated. Tables 14 and 15 compile these 
results for a variety of Reynolds numbers and locations along the cylinder 
boundary. 


Table 14. Boundary layer thickness for Oseen flow. 


12 0/tt 

R = 2.0 

R = 3.0 

R » 4.0 

R = 5.0 

A(9)/r'Y 

12 

0. 357 

0. 236 

0.169 

0.125 

0. 798/r1*13‘+ 

11 

0. 365 

0.243 

0.175 

0.131 

0. 804/r1*113 

10 

0.390 

0. 264 

0.194 

0.148 

0. 823/ R^* 0 5^ 

9 

0. 436 

0. 301 

0.227 

0.178 

O. 865/r0*9 72 

8 

0. 509 

0. 361 

0. 280 

0. 227 

0. 941/r0*880 

7 

0. 619 

0. 451 

0. 360 

0. 300 

1. 071/rO*788 

6 

0. 786 

0. 589 

0.482 

0. 412 

1. 28O/r8*7 0*+ 

5 

1.048 

0. 807 

0. 674 

0. 589 

1. 618/r0*6 30 

4 

1.484 

1.169 

0.996 

0. 885 

2.188/ r0*565 

3 

2.278 

1. 833 

1.589 

1.431 

3. 226/r0*5 08 

2 

3. 990 

3. 276 

2. 883 

2. 629 

5. 449/R8*‘+57 

1 

9. 414 

7. 876 

7. 028 

6. 478 

IZv 438/R°*'+03 
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Table 15. Boundary layer thickness for Burgers flow. 


12 0/tt 

R » 2.0 

R - 2.83 

R " 3.46 

R ■ 4.0 

A(0)/rY 

12 

0.574 

0.482 

0.441 

0.410 

0.800/rO‘‘*32 

11 

0.581 

0.488 

0.445 

0.415 

0.811/R°*'+3'+ 

10 

0.602 

0.505 

0.460 

0.428 

0.843/r0*‘*89 

9 

0.639 

0.535 

0.487 

0.453 

0.900/rO’‘^3^ 

8 

0.699 

0.582 

0.529 

0.491 

0.991/R°*®°^ 

7 

0.789 

0.654 

0.593 

0.549 

1.130/R°‘520 

6 

0.927 

0.765 

0.690 

0.638 

1.342/R°-S37 

5 

1.144 

0.938 

0.844 

0.778 

1.679/R°*5^® 

4 

1.510 

1.231 

1.102 

1.014 

2.244/R°'^'^^ 

3 

2.190 

1.780 

1 . 588 

1.457 

3.285/R°*^®^ 

2 

3.694 

3.007 

2,681 

2.459 

5.541/R°-586 

1 

8.579 

7.050 

6.316 

5.814 

12.639/R°‘^®° 


In both cases the thickness of the boundary layer increases as one 

moves around the cylinder from the forward stagnation point (0 * ir) . 

Increasing the Reynolds number serves to compress the boundary layer. The 

last column in each table is a least- squares fit of the data given in each 

row to an expression of the form A(0)/R^. The function A(9) is similar 

in both flows. The value of y depends on 6 in the case of Oseen flow* 

but appears to hover about the constant — in Burgers flow in the range 

2 

JL < 9 < TT. This suggests that boundary layer thickness in Burgers flow is 
2 

roughly proportional to 

Analytic studies of the boundary layer on a serai-infinite flat plate 
using the full nonlinear equations show that the thickness is proportional 
to Our work, although not conclusive, suggests a similar result for 

a circular cylinder using a linear model. The apparent agreement between 


2,6 



i/ 




this prediction of Burgers flow and that of nonlinear analysis verifies 
agvfin its superiority over Oaeen flow in describing fluid behavior near the 
cylinder. 


k 
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APPENDIX I 


The standard form of Mathieu's equation is, from McLachlan (1964), 

y" 4- (o - 2q cos 26 ) y " 0 (1*1) 

where q is termed the modulus and a is the eigenvalue. In this 
discussion we are only interested in solutions of (1. 1) which are odd in 0#. 
In many applications the eigenvalue is determined from the condition that 
y(0 ) be periodic with period 2ir. Thus if q “ 0, the eigenvalue must be 
the square of an integer (a^ » ti? ) and the corresponding odd 
eigenfunction is 

y^(6 ) “ sin n0 (1*2) 

If Iql is small but nonzero, the eigenvalue and eigenfunction can be 
expanded in series of the form 

°n “ 

k“l 

00 

yj^(9»q) = sin n0 + ^ 

k“l 

The coefficients and ) sre determined by substituting these 

expansions into (I.l), equating like powers of q, and requiring that 

S^k(6 ) bs in 6 and periodic with period 2tt. The expression for 

y (0 ) can be rewritten in the fom 
n 

(1.5) 

A normalized tlathieu function 


OO 

=■ I D^"\q) sin k0 

il .l ^ 


k-1 


where each Di|^'^\q) is a power series in q. 
is defined as follows: 
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00 

se^(e,q) - y^(0,q)/[ I {Dj^"\q)}2]l/2 

k“l 

00 

= I sin k0 

k“l 

00 

where I [si^\q)]2 = i, 
k=l 

'The standard form of the modified l-lathieu equation is 

Y"(z) “ [<^ 2q cosh 2z] Y (z) » 0 
n n n 

IJnen q < 0, the eigenfunctions which decay exponentially as z 
denoted by 

Y^(z) = Gek^(z,q) 


From McLachlan (1964), we have 


Gek (z,q) 
n 


c K ( Iq e^) 

n n^ ^ '' 

c(i.ir) exp [-JLz - e^] 

ri 2 2 


( 1 . 6 ) 


(1.7) 

“ are 

( 1 . 8 ) 


(1.9) 


where q < 0. 
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APPENDIX II 


The vorticity equation for Oseen flow past a circular cylinder can be 
solved using the method outlined in section 3. The appropriate coordinate 
systera (?,n) is defined by 


5 

X “ e cos ri, 


y “ e sin n 


(II. 1) 


and the metric coefficient is 


ot(C,n) “ 


The vorticity is given by 


= F(?,n) exp [— R e^ cos n] 

9 


where F(5ih) satisfies the equation 


F = 0 

9^2 9ti2 4 


(II. 2) 


(II. 3) 


(II. 4) 


Solutions to this equation are of the form 

F (5,n) = E (5) sin nn 
n n 

» 

where En(5) satisfies the equation 

E" -(n2 + 1 r2 E = 0 

n , n 


From Gradshteyn and 


E (O = A I 
n n n 


Ryzhik (1965), we have 
(i R e') . (i R e«) 


(II. 5) 


(II. 6) 


(II. 7) 
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since the vorticlty must decay exponentially as ^ the coefficients 

Ajy are all zero and the vorticlty is given by 

00 

,n ) “ -exp [i. R e^ cos n] I B K (i R e^) sin nh (II. 8) 

2 n-1 " 2 

In terms of polar coordinates (r,Q), the asymptotic behavior of the Oseen 
vorticlty is 

m (r,9 ) w(9 ) r~^^^ exp [ •" il Rr (1 - cos 0 )] (II. 9) 

2 

where w(0) = 0, w’ (0) * 0 and w(9 ) has period 2 ti. 
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APPENDIX III 

The expansion of 4' in (7. 8) contains the expression 

lim (C »h ) w (? ,n ) sin kn dn (III. 1) 

5-Ho 0 

which can be computed once the asymptotic behavior of and w is known. 

From (II. 2), (3.24) and (3.30), we have 

a2 (? ,T1 ) ~ e^^ (III. 2) 

We also know that 

00 

w(c ,n) “ - exp [AR cosh ? cos D ] W C3ek (? , - i A^ R^} 

n=l " 4 

• sejh , - i A^ R^) (III. 3) 

4 


From (1.9) the asymptotic form of the vorticity is given by 


/ir \l/2 T 1 r , 

w (C >n ) *^ •“ ( — 1 exp [ ” — ? - — AR e (1 - cos n )] 
\AR/ 2 2 


• I c W se (n , - i a2r2) 
n n o'* ' . 

n=l 4 


Thus the integral in (III. 1) simplifies to 


( 3 \ 1 / 2 

exp [ - i (AR e^ - 35 ) ] 
R / 2 

00 

* \ G W exp [i AR e^ cos n] sin kn 

n“l o 2 

• se (n , - i A^R? ) dri 

" 4 


(III. 4) 


(III. 5) 
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From McLachlan (1964), we have 

00 

se (h, - “ A^r2) “ 5] S ^ sin (2m + e )n 

4 m-0 ” 

where e * 1 if n is odd, 
n 

* 2 if n is even. 

Therefore, we have 


se (ri, “ — A^R^) sin kri “ — I S [cos (2m + e - k)Ti 
^4 2 m»0 "* " 


- cos (2m + e + k)ri] 
n 


(III. 6) 


(III. 7) 


The substitution of (III. 7) into (III. 5) yields integrals of the form 

exp [u cos 0] cos n0 d0 » ir I^^Cu) (ill. 8) 

o 


where is the modified Bessel function of the first kind. After 

evaluating the integral in (III. 5) and replacing the Bessel functions by 
their asymptotic forms, we have 


0) sin kndn ~ 
o 


— k I c I (2m + e ) + 0(e"^) (III. 9) 

r 2 n=l ^ " m=0 " ® 


Taking the limit in (lll.l) and substituting back into (7.8), we obtain 


^00 ^ 
ijj(?*,ri') ~ A e^ sin n' - — C_ ^ — sin kri' + 0(e ^ ) 

TT “ k=l k 


(III. 10) 


where 


C = 
D 


_ 4ir 


y c W y (2m + e ) S 
n n n 1 

R'^ n=l m=0 


(n) 


(III. 11) 
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But 


00 

— sin kn' - i(±TT - n’) (III. 12) 

k"l k 2 

where the plus sign is chosen when 0 < H ’ ^ 3^nd the minus sign when 

-ir < n ' < 0. Thus the asymptotic form of the stream function for Burgers 
flow past an elliptical cylinder is given by (7.9). 

This work was supported by the National Aeronautics and Space Adminis” 
tration under contract no. NAGl-197. 
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Figure 3 
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Convective velocity field at the rear stagnation point of the 
cylinder in Oseen flow. 
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Figure 4 
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Convective velocity field in Burgers flow. 
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Figure 7. Drag coefficient vs. Reynolds number for Burgers flow past a 

variety of geometries. , flat plate perpendicular to 

flow; - ~ , ellipse (aspect ratio 3:1) perpendicular to 

flow; ; , circle; , ellipse (aspect 

ratio 3:1) parallel to flow; — — , flat plate 

parallel to flow; 0 , numerical solution (Dennis and 
Dunwoody) of the full nonlinear equation for flow past a flat 
plate. I 

ft' 

ij 
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